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Abstract Papilio machaon is distributed in the subarctic and temperate zones of the northern hemisphere, including 
the Eurasian and North American Continents. It is also distributed in the Japanese Islands and Sakhalin, and is 
classified as the subspecies hippocrates. In order to elucidate the phylogenetic relationship between the Japanese P. 
machaon population and Continental populations, a molecular phylogenetic analysis was performed with 
mitochondrial NDS DNA sequences using P. machaon specimens collected from various areas in the Japanese Islands 
and foreign countries and other species included in the subgenus Papilio. We found that the Japanese P. machaon 
population (the subgenus hippocrates) was genetically distinct from the Eurasian and North American populations. 
The Japanese population diverged earlier than other Continental P. machaon populations in the subgenus Papilio, 
which indicates that the Japanese population would be isolated in the Islands since their geographical establishment. 
These results imply that the Japanese population of other butterfly species may also be distinct from the Continental 


populations at the molecular level even though morphological similarities exist between the populations. 
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The genus Papilio Linnaeus, 1758 consists of various 
morphologically distinct groups, some of which have been 
further classified into subgenera such as Papilio, Euchenor, 
Menelaides, and Achillides (Igarashi, 1979). The subgenus 
Papilio includes approximately 10 species that are distributed in 
the northern hemisphere (Fujioka et al., 1997). Papilio 
machaon, P. hospiton, P. saharae, and P. sikkimensis are 
distributed in the Eurasian Continent and its affiliated areas 
(Fujioka et al., 1997), while various species such as P. machaon, 
P. polyxenes, P. zelicaon, P. indra, P. brevicauda, P. hudsonianus, 
P. oregonius, and P. bairdii are distributed in the North American 
Continent (Scott, 1986; Fujioka et al., 1997). The classification 
of these species is complex because some hybridize in their 
common habitats (Dupuis and Sperling, 2015; 2016), and some 
of them are regarded as subspecies by some authors (Tyler, 
1994 ; Fujioka et al., 1997). Papilio machaon is widely 
distributed in the subarctic and temperate zones from Europe to 
northern America through Russia, Mongolia, China, Alaska, and 
Canada, and many subspecies have been recorded in various 
regions. In east Asian islands including Japan, P. machaon is 
distributed in Sakhalin, Hokkaido, Honshu, Shikoku, Kyushu, 
and their affiliated small islands (Fujioka et al., 1997), and this 
population has been classified as the subspecies hippocrates C. 
et R. Felder, 1864. 
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Papilio machaon is one of the common papilionid butterflies in 
the Japanese Islands, including Sakhalin, and is distributed in 
wide areas, except for the southern Ryukyu islands. The 
difference in latitude between the northernmost and 
southernmost habitats, Sakhalin and Yakushima Islands, 
respectively, is approximately 24 degrees. According to the 
records of the Japan Meteorological Agency (2017), the 
difference in the average annual temperature between these 
islands is approximately 17°C. Larval host plants of P. machaon 
are various local species of Umbelliferae. In the middle area of 
Honshu, their larvae were found at the seashore as well as in 
mountains with an altitude of 2,000 m. These findings indicate 
that P. machaon is adapted to highly variable climate conditions 
in the Japanese Islands. 


Previous studies elucidated the molecular phylogenetic 
relationships among species of the subgenus Papilio (Spering 
and Harrison, 1994; Aubert et al., 1999; Caterino and Sperling, 
1999 ; Reed and Sperling, 1999 ; Zakharov et al., 2004); 
however, the Japanese population was not examined in detail. 
Many types of endemic insects including butterflies such as 
Luehdorfia japonica, Lethe sicelis, and Neope goschkevitschii 
inhabit the Japanese Islands (Kawazoe and Wakabayashi, 1976; 
Shir6zu, 2006). In some butterfly species inhabiting the 
Japanese Islands and Eurasian Continent, e.g., Parnassius 
citrinarius, Artogeia napi, and Erebia niphonica, the Japanese 
population is largely divergent from the Continental populations 
at the molecular level (Yagi et al., 2001 ; Fujii et al., 2001 ; 
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Nakatani et al., 2007 ; Shinkawa and Nonaka, 2010). These 
differences in molecular phylogeny may reflect the geographical 
isolation of the Japanese Islands since their formation (Tojo and 
Ito, 2015). The repeated connection and separation of the Islands 
at the eastern edge of the Eurasian Continent during the late 
Pliocene and Pleistocene resulted in the unique characteristics of 
butterfly fauna in the Japanese Islands (Tojo and Ito, 2015). 


Papilio machaon is a widely distributed species in the northern 
hemisphere, whereas the Japanese Islands population has unique 
characteristics in terms of climate adaptation. Genetic differences 
in the Japanese population from those in other areas are of great 
interest to entomologists. In the present study, P. machaon 
individuals were collected at various areas in the Japanese 
Islands including Sakhalin, and the nucleotide sequences of their 
mitochondrial NADH dehydrogenase subunit 5 (ND5) genes 
were elucidated. Molecular phylogenetic trees of P. machaon 
individuals from the Japanese Islands and foreign areas were 
constructed with other species of the subgenus Papilio in order 
to clarify phylogenetic relationships with the Japanese 
population in the subgenus Papilio. 


Materials and methods 


Sample collection 


Specimens of P. machaon hippocrates were collected in various 
areas of Hokkaido, Honshu, Shikoku, Kyushu, and their 
affiliated small islands in Japan between 2006 and 2009. Foreign 
specimens were collected between 1996 and 2003. Adult bodies 
without wings and larvae were preserved in 99.5% ethanol at 
-20°C. Foreign samples that we had previously described (Yagi 
and Fujioka, 2005) were used for molecular phylogenetic 
analyses in the present study. Specimens from Japan and foreign 
countries were kindly provided by many entomologists. Species, 
sampling locations, and database accession numbers are shown 
in Table 1. 


DNA extraction, ND5 gene amplification, and sequencing 


Adult legs or larval muscles were shredded by scissors in a test 
tube, and mixed with 500 ul lysis solution containing 0.1% SDS 
and 0.5 mg/ml proteinase K. The test tube was then kept on a 
rotator for 30-120 min at 55°C. The digested samples were 
mixed with phenol once, phenol/chloroform (1 : 1) twice, and 
chloroform/isoamylalcohol (24:1) once. DNA in the water 
fraction was precipitated with ethanol, washed twice with 75% 
ethanol, and dried. Dried DNA was resuspended in 50 ul Tris- 
EDTA buffer (pH 7.5) and stored at -20°C. 


A part of the NDS gene of mitochondrial DNA was amplified 
with a polymerase chain reaction (PCR) using the following 
primers: V 1 (5’-CCTGTTTCTGCTTTAGTTCA-3)), K3 (5’ 
-TAKCTTCAATATTAYRCTCT-3’) (Yagi et al., 1999), and 


KIAGEHA (5'--TAGGACAAAGTTTATTAAAG- 3’). The PCR 
reaction mixture contained 25 ul buffer, 10 ul dNTPs, 1.5 ul V1 
primer (10 uM), 1.5 ul K3 primer (10 uM), 10 wl sterilized 
water, 1 ul butterfly DNA (~4 ng), and 1 ul DNA polymerase 
(KOD FX, TOYOBO, Osaka, Japan). PCR reactions were 
performed using the following amplification profile: 94°C for 2 
min, followed by 30-40 cycles at 98°C for 10 sec, at 45-50°C for 
30 sec, at 68°C for 1 min, and final extension at 68°C for 5 min. 
The PCR product was inspected by 1% agarose gel 
electrophoresis, and purified using the Sephacryl S-300 High 
Resolution column (GE Healthcare, Little Chalfont, UK). 


A DNA sequencing reaction of the PCR product was performed 
with the same primers for PCR and the BigDye Terminator 
Cycle Sequencing Kit (Applied Biosystems, Waltham, MA). 
DNA sequences were obtained with the ABI 3100 Genetic 
Analyzer (Applied Biosystems, Waltham, MA). 


DNA primers were synthesized by Sigma-Aldrich Japan, Tokyo, 
and other chemicals were purchased from Nacalai Tesque, 
Kyoto; TAKARA BIO, Kusatsu; and Sigma-Aldrich Japan, 
Tokyo. 


Phylogenetic analysis 


ND5 gene nucleotide sequences (717 bp) were aligned using 
Clustal W and fragment comparisons were conducted using 
BioEdit version 7.0.9.1 (Hall, 1999). A phylogenetic analysis of 
P. machaon samples obtained at various locations in the Japanese 
Islands and other countries was conducted using the Maximum 
likelihood (ML) method and Neighbor-joining (NJ) method 
with MEGAS software (Tamura et al., 2011). In the analysis, 
evolutionary history was inferred using the Tamura-Nei model 
(Tamura and Nei, 1993) with a uniform base substitution rate 
and bootstrap of 10,000 replicates. The codon positions included 
were Ist+ 2nd+ 3rd+Non-coding. The base substitution was 
confirmed to be unsaturated by the function included in MEGAS. 
The ML tree was drawn with the highest log likelihood values 
and a scale of the branch length proportional to the number of 
substitutions per site. The NJ tree was drawn with the values of 
bootstraps and branch lengths higher than 0.003. 


Results and Discussion 


The ML and NJ molecular phylogenetic trees of the subgenus 
Papilio were shown in Figures 1 and 2 , respectively. When P. 
xuthus and P. benguetana were settled as outgroup species, P. 
indra diverged first; the pair of P. polyxenes and P. zelicaon 
diverged second; the pair of P. brevicauda and P. hudsonianus 
diverged third; and the P. machaon group finally diverged. 
Among P. machaon, it is important to note that the Japanese 
population (clade V: subspecies hippocrates) diverged first. The 
divergence order of other P. machaon populations currently 
remains unclear due to low bootstrap values; however, P. 
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Table 1. Papilio samples used in the phylogenetic analysis. 


Species 


Accession 


number 


Location collected 


123 





P. machaon 


LC061680 


Japan: Hokkaido: Esan, Honshu: Omonogawa, Gojome, Oga, Kisakata, 
Sabusawashima, Katsurashima, Shichigahama, Yuza, Sagae, Mito, Tainai, Ito, 
Shimoda, Tatsunokuchi, Hakusan, Mikata, Joyo, Awajishima, Maniwa, 
Yamaguchi, Shikoku: Miyoshi, Marugame, Manno, Kyushu: Iizuka, Kanzaki, 
Aso, Tsushima; Korea: Jeju-do 





P. machaon 


LC061681 


Japan: Honshu: Oma, Oga, Kisakata, Sagae, Mito, Joyo, Maniwa, Shikoku: 
Miyoshi, Manno, Kyushu: Iizuka, Aso 





P. machaon 


LC061665 


Japan: Hokkaido: Wakkanai, Kushiro, Furano, Esan, Rishiri, Rebun; Russia: 
Sakhalin 





P. machaon 


LC061725 


Japan: Hokkaido: Wakkanai, Furano, Esan, Honshu: Oma; Russia: Sakhalin, 
Magadan 





P. machaon 


LC061717 


Russia: Sakhalin 





P. machaon 


LC061684 


Japan: Honshu: Shimoda, Hakusan, Joyo 





P. machaon 


LC061706 


Japan: Honshu: Kisakata, Shimoda, Yamaguchi 





P. machaon 


LC061720 


Japan: Honshu: Oma, Sabusawashima, Sagae, Mito, Tainai, Hakusan 





P. machaon 


LC061736 


Russia: Magadan, Chita, Kamchatka, Tuva, East Sayan; Tajikistan: Pamir Mts; 
Mongolia: Terelj, Tsenkher; China: Gansu; Macedonia; USA: California, 
Oregon: Spensbridge; Canada: Yukon 





P. machaon 


LC061737 


Russia: Kamchatka 





P. machaon 


LC274577 


USA: California: San Bernardino 





P. machaon 


LC274578 


Russia: East Sayan 





P. machaon 


LC061738 


Russia: Chita 





P. machaon 


LC061739 


Russia: Primorsky Krai 





P. machaon 


LC274579 


Russia: Altai 





P. machaon 


LC274580 


Canada: Yukon 





P. machaon 


LC061740 


China: Yunnan: Mt Guangya 





P. machaon 


LC061730 


Korea: Gyeonggi-do: Yonjongdo 





P. machaon 


LC061733 


China: Yunnan 





P. machaon 


LC274581 


Russia: Altai 





P. machaon 


LC274582 


Tajikistan: Khodza-Mumin 





P. machaon 


LC274583 


Russia: Orenburg 





P. machaon 


LC061734 


Russia: Primorsky Krai 





P. machaon 


LC061735 


China: Beijing 





P. saharae 


LC274584 


Morocco: Tezenakht 





P. saharae 


LC274585 


Morocco: Tezenakht 





P. sikkimensis 


LC061732 


China: Tibet 





P. sikkimensis 


LC274586 


China: Tibet 





P. sikkimensis 


LC274587 


China: Tibet 





P. brevicauda 


LC274588 


Canada: Newfoundland 





P. hudsonianus 


LC274589 


Canada: Hudson Bay 





P. polyxenes 


LC061731 


Canada: Ontario: Komoka, Manitoba: Duck Mountain 





P zelicaon 


LC274590 


USA: Idaho 





P. indra 


LC274591 


USA: Colorado: Sweet Pass 





P. benguetana 


LC274592 


Philippines: Benguet 





P. xuthus 








LC274593 





Korea: Gyeonggi-do: Yonjongdo; Japan: Honshu: Kyoto 


*indicates the clade of P. machaon phylogenetic trees (Fig. 1 and 2). 
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Fig. 1. ND5 phylogenetic tree of the subgenus Papilio including P. machaon collected at various 
locations in the Japanese Islands and other countries, constructed with the Maximum likelihood 
(ML) method using MEGAS software. The time scale is based on the tentative evolution rate 
estimated from Papilio binor ND5 (0.81% substitutions/site per m.y.) (Osozawa et al., 2013). 
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Fig. 2. ND5 phylogenetic tree of the subgenus Papilio including P. machaon collected at various 


locations in the Japanese Islands and other countries, constructed with the Neighbor-joining 
(NJ) method using MEGAS software. 
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machaon was clearly divided into 5 clades. Clade I comprised 
25 individuals of the Eurasian and North American Continents 
with 9 haplotypes. Clade H included 9 individuals of Eurasia 
with 7 haplotypes. Clade III contained 6 individuals of P. 
saharae from northern Africa. Clade IV comprised 3 individuals 
of P. sikkimensis from Tibet. Clade V comprised 312 individuals 
of P. machaon from the Japanese Islands, Sakhalin, and Jeju-do 
Island (Korea) with 8 haplotypes. A similar phylogenetic 
topology was obtained by the ML and NJ methods (Fig. 1 and 
2). In the Eurasian Continent, Clade I and II individuals were 
distributed sympatrically. As reported by Fujioka et al. (1997), 
P. saharae and P. sikkimensis are different species from P. 
machaon, but are included in P. machaon in these phylogenetic 
trees (Fig. 1 and 2). 


The first divergence of P. indra among the subgenus Papilio in 
our NDS trees is consistent with previous findings in which the 
ML and MP (maximum parsimony) phylogenetic trees were 
constructed using mitochondrial ND/ and COI/COI genes 
(Aubert et al., 1999 ; Caterino and Sperling, 1999 ; Reed and 
Sperling, 1999), and the pairing of P. zelicaon with P polyxenes 
in our ML and NJ trees was also reported in previous studies 
(Caterino and Sperling, 1999 ; Reed and Sperling, 1999); 
however, the divergence order of P. indra, the P. zelicaon and P. 
polyxenes pair, P. hospiton, and P. machaon was controversial. 
In recent studies by Sperling and his colleagues (Zakharov ef 
al., 2004 ; Dupuis and Sperling, 2015), they constructed more 
concrete ML and Bayesian phylogenetic trees using the COJ/ 
COI genes with P. xuthus as an outgroup species, which showed 
that divergence occurred in the order of P. indra, the P. zelicaon 
and P. polyxenes pair, P. hospiton, and finally the P. machaon 
species group. Our present results support this divergence order. 
Papilio hospiton is not included in the present study, but 
discussed in more detail below. 


All samples of North American P. machaon examined were 
involved in Clade I (Fig. 1 and 2). In the North American 
Continent, P. machaon has been classified into many species/ 
subspecies, e.g., P. m. aliaska, P. bairdii, P. brevicauda, P. m. 
dodi, P. hudsonianus, P. joanae, P. m. kahli, P. m. oregonius, and 
P. m. pikei, among which genetic introgression among some 
species/subspecies has been reported (Dupuis and Sperling, 
2015; 2016). Based on the results of the present study, we were 
unable to elucidate phylogenetic relationships among North 
American P. machaon because of the small number of analyzed 
samples. 


Regarding P. m. hippocrates, Sperling and his colleagues used 
only a few samples from Japan in their phylogenetic analysis, 
and showed that it diverged earlier than other P. machaon 
subspecies from France and Washington, USA (Sperling, 1993; 
Zakharov et al., 2004). Our present (clade V) and early pilot 
phylogenetic studies (Yagi and Fujioka, 2005) using a large 
number of P. m. hippocrates samples from various areas of the 


Japanese Islands and Sakhalin supported the findings of 
Zakharov et al. (2004). P m. hippocrates was clearly distinct 
from other P. machaon subspecies at the molecular level, 
suggesting vicariant speciation due to its isolation in the 
Japanese Islands for a few million years. It is currently 
inconclusive whether the Japanese population is an independent 
species from Eurasian and North American P. machaon due to 
the difficulties associated with massive hybridization 
experiments; however, Ae (1971; 1988) reported an unbalanced 
F; sex ratio and low F; fertility after the subspecies-crossing of 
Japanese hippocrates with German gothica and British 
britanicus. Ae (1988) found that the eggs from the F, 
(hippocrates x britanicus) mutual cross (F, eggs) and F, male 
back-cross with female hippocrates rarely hatched. It is 
important to note that Remington (1959) raised the subspecies 
hippocrates to species because of differences in morphological 
characteristics and chromosome numbers from European P. 
machaon and a marked deficiency in females in the broods of F, 
hybrids with European P. machaon, The major chromosome 
numbers of Japanese and European P. machaon were previously 
reported to be 62 and 60, respectively (Maeki, 1953 ; 1976). 
Phylogenetic analyses performed by Sperling and our groups in 
addition to previous findings clearly indicated that the Japanese 
P. machaon population (hippocrates) shows large genetic 
differences from the other populations. 


Papilio hospiton is endemic to the Mediterranean Sardinia and 
Corsica Islands, which are located near the Italian Peninsula. 
Although P. hospiton was not used in the present analysis, 
previous studies showed that it diverged earlier than all P. 
machaon subspecies in phylogenetic trees of the subgenus 
Papilio (Sperling, 1993 ; Zakharov et al., 2004 ; Dupuis and 
Sperling, 2016). Papilio hospiton was found to be an 
independent species from P. machaon based on hybridization 
and the findings of molecular studies; however, the low 
frequency of genetic introgression has now been observed 
between the two species in the field (Cianchi et al., 2003). 
Papilio hospiton may have been isolated in the Sardinia and 
Corsica Islands for a long time without genetic exchange with 
Continental P. machaon. A similar relationship between Japanese 
(hippocrates) and Continental P. machaon has been suggested. 


The present Papilio phylogenetic trees using mitochondrial ND5 
gene may not be consistent entirely with speciation history of 
this group because genetic introgression occurs among the North 
American species/subspecies group (Dupuis and Sperling, 2015; 
2016). Although P. hudsonianus (V) is generally recognized as 
a subspecies of P machaon (Scott, 1986), it clustered with P. 
brevicauda (U) (Fig. 1 and 2), which may be due to the 
mitochondrial DNA introgression. Future work with nuclear 
genes, either gene sequences, microsatellites, or genome-wide 
SNPs, will help to elucidate the species phylogeny in this group 
(Dupuis and Sperling 2016). 
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In conclusion, we found that the Japanese P. machaon population 
classified as the subspecies hippocrates is genetically distinct 
from the Eurasian and North American populations. Japanese P. 
m. hippocrates diverged earlier than other Continental P. 
machaon subspecies within the subgenus Papilio, which 
indicates that the Japanese population was isolated in the Islands 
since their geographical establishment. Molecular variations in P. 
machaon within the Japanese Islands need to be analyzed further 
in order to elucidate phylogeography in more detail. The results 
of the present study imply that the Japanese populations of other 
species are also distinct from the Continental populations at the 
molecular level even though their morphological characteristics 
are similar between populations. 
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FTF I YT REB LUKT AU AKO WEY 
B ABACK l DATA. FTH RNB k 
OFANY VICS DAL, Zn Ò lE hippocrates & LT 
KMESNA. FPTFANONASI S (FAV eG) KH 
CKHRMORRMRe HDMF ARO, AA 
MOET TNBEL OSE OMIR 5 DNA ZABEL, 3 F 
IY EY DNA ®D NDS tte TO TBH AR ALY & PEL, 
TTAR eH Te. CORR, HAVE RAILL—7Y 
TREN B k UIT 2 Y ASAD 5 HRY IC HATE IC 
SNACEBHPoK. FT TNODLARRBICBWT, 
AAW BEAL —-FYV PREM LILT AY HENK Y 
JNDPMCWR. COLE, FT PVA, KEP 54 
PUAAARIBICBV TC, RMA & yii l WAH 
SHRTOKCESRBLTWS. HAVIN Kk 
& DE CHEDDW Cv CT d HRN ICHEN TVS FEDS, fC 
bFTET SDDS LANAY. 
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